ASYMPTOTIC BEHAVIOR FOR A NONLOCAL DIFFUSION EQUATION 
WITH ABSORPTION AND NONINTEGRABLE INITIAL DATA. THE 

SUPERCRITICAL CASE 



JOANA TERRA AND NOEMI WOLANSKI 



Abstract. In this paper we study the asymptotic behavior as time goes to infinity of the 
solution to a nonlocal diffusion equation with absorption modeled by a powerlike reaction — u p , 
p > 1 and set in M. N . We consider a bounded, nonnegative initial datum uo that behaves like a 
negative power at infinity. That is, \x\ a uo(x) — > A > as \x\ — > oo with < a < N. We prove 
that, in the supercritical case p > 1 + 2/a, the solution behaves asymptotically as that of the 
heat equation -with diffusivity a related to the nonlocal operator- with the same initial datum. 



1. Introduction 

In this paper we study the asymptotic behavior as time goes to infinity of the solution to a 
nonlocal diffusion equation in R N with an absorption term. We are interested in the case in 
which the initial datum is nonintegrable. More precisely, we consider a bounded datum uq(x) 
such that |x| Q uo(x) — > A > as \x\ — > oo with < a < N. The equation under consideration is 
the following 

(1.1) u t (x,t)=J J(x - y)(u(y,t) - u(x,t)) dy - u p (x,t) 

where J € C^ C (K N ), radially symmetric with J > and f J = 1. 

The term f J(x — y)(u(y,t) — u(x,t)) dy in (jl.ip represents diffusion. In fact, let u(x,t) be 
the density of a population at the point x at time t. The kernel J(x — y) may be seen as the 
probability distribution density of jumping from site y to site x. Thus, f J(x — y)u(y, t) dy 
represents the rate at which the population is arriving at the site x from all over space. By the 
symmetry of the kernel J, J J(x — y)u(x, t) dy represents then the rate at which it is leaving the 
point x. 

The absorption term —u p (x,t) in the equation represents a rate of consumption due to an 
internal reaction. 

This diffusion operator has been used to model several nonlocal diffusion processes in the last 
few years. See for instance [TJ [2j [31 HI [19]. In particular, nonlocal diffusions are of interest 
in biological and biomedical problems. Recently, this kind of nonlocal operators have also been 
used for image enhancement [10] , 



It can be seen that -when properly rescaled- the nonlocal operator converges to aA, where A 
is the Laplace operator and a is a constant that depends only on J and the dimension N. (See, 
for instance, [3[H]). 
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Another close relation to the heat operator with diffusivity o was found in [51 [13] where it 
was proven that, for bounded and integrable initial data, the asymptotic behavior as t tends to 
infinity of the solution ul to the equation without absorption 

(1.2) u t (x,t) = j J(x- y)(u(y,t) - u(x,t)) dy = Lu, 

is the same as that of the solution of the heat equation with diffusivity a and the same initial 
condition. Namely, 

(1.3) t N/2 \\u(x,t) - MU a {x,t)\\ Lao{lxl < KVi) ^0 ast^oo VK > 0, 

with U a the fundamental solution of the heat equation with diffusivity a and M = j uq the 
initial mass. The assumption on the initial condition is that uq £ L°° n L . 

In a recent paper [T7], the nonlocal equation with absorption (jl.ip was considered in the case 
u(x, 0) G L°° n L 1 and p > 1 + 2/N. The authors prove that u also satisfies (|1.3p . Some results 
for p < 1 + 2/N can also be found in that paper. The critical case p = 1 + 2/N was left open. 

On the other hand, the asymptotic behavior of the solutions to the heat equation and of the 
heat equation with absorption are very well known for some nonintegrable initial data with a 
precise power like behavior at infinity. In fact, Kamin and Peletier in [JJj study the asymptotic 
behavior of the solution to 

u t = aAu-u p in R N x (0,oo), 

(1-4) ,r 

u(x,0)=u (x) inR N , 
with no £ L°° and uq > such that 

(1.5) |x| a ?/o(x) — > A > as \x\ — > oo, 
with < a < N. In fact, the authors prove that 

(1.6) t a / 2 \\u(x,t)-U a)A (x,t)\\ LOC(l , 
where U aj A is the solution to 

U t = aAU - 5 P:a U p in R N x (0, oo), 



(1-7) „, A 



J p,a'- 

U(x,0) = -j- in 



Here 5 Pta = in the supercritical case p > 1 + 2/ a and 5 Pt0l = 1 in the critical case p = 1 + 2 /a. 

Then, in [15] . Kamin and Ughi study the case where a = N in (|1.5p and p > 1 + 2/iV. They 
prove that 

(1.8) t N > 2 ^-C AN U a (x,t) 



-+ as t ->• oo VA' > 0, 



with C^at a constant that depends only on ^4 and A and U a the fundamental solution of the 
heat operator with diffusivity a. 

Later on, in [12] , Herraiz fully analyzes problem (jl.4p and presents a complete characterization 
including the asymptotic behavior outside the parabolas \x\ < Ky/t. 



It is the purpose of this paper to prove that, for initial data satisfying f 1 1 . 5 j) with < a < N 
and p > 1 + 2/ a, the asymptotic behavior of the solution to the nonlocal diffusion problem with 
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absorption (jl.ip is the same as the asymptotic behavior of solutions to the heat equation with 
absorption (jl.4p . That is, that (|1.6p holds in the case a < N and (jl.8p holds in the case a = N 
with u the solution of (jl.ip . 

The proofs of these results rely strongly on estimates for the fundamental solution of the 
equation without absorption (|l,2p . 

In [3] the authors observed that the fundamental solution of the operator in (jl.2p can be 
written as 

U(x,t) = e- l 8 + W(x,t), 

where 5 is the Dirac measure and W(x, t) is a smooth function. 

Then, in [T2j the authors established L°° estimates for W and proved that 

t N/2 \\W(x,t) - J7 a (x,i)|| i0 o( K JV) ->• as t ->■ oo, 

where U a is the fundamental solution of the heat equation with diffusivity a. 

This estimate allowed them to prove that the solution of the homogeneous nonlocal diffusion 
equation (jl.2p with bounded and integrable initial data behaves -for t tending to infinity- as the 
solution of the heat equation with diffusivity o, removing in this way the hypothesis that the 
Fourier transform of uq be in L 1 as was assumed in [5]. 

In the present paper, as we are considering nonintegrable initial data, these L°° estimates 
cannot be used. So, the first step is to get estimates in L 9 — norms for 1 < q < oo. 



As in [T7], the philosophy of the study of the asymptotic behavior for the problem with 
absorption in the supercritical case is that, as time tends to infinity, the absorption term goes 
to so fast that the behavior is the same as that of the equation without absorption. Thus, 
what we need to prove is essentially that, as time tends to infinity, the difference between the 
solutions to both problems goes to faster than each term separately. 

In order to prove this result it is necessary to compare the solution u of (jl.ip with the solution 
ul of the equation without absorption (|1.2p that coincides with u at a large time to- Once we 
prove that this difference is small the result follows if we know that the asymptotic behavior of 
ul is the same as that of the heat equation with diffusivity a. To this end, we first analyze the 
asymptotic behavior for the nonlocal equation without absorption (|1.2p with an initial datum 
satisfying (|1.5p . In order to finish the argument, it is essential to know, and we prove it in this 
paper, that the behavior in space of the solution of our problem (|1.1|) at any time to is the same 
as that of the initial datum no- That is, that u(-,to) satisfies (|1.5p at any positive time to- 

The techniques of this paper do not work in the critical case p = 1 + 2/a. In fact, for such p, 
the asymptotic behavior is that of the heat equation with absorption since the diffusion and the 
absorption terms are of the same order. We have to proceed in a completely different way by 
going back to the techniques of the original paper for the heat equation p3], and we do so in |18j . 
The method in [18] also works in the supercritical case considered in this paper. Nevertheless, 
the ideas in the present paper are simpler and rather easy to implement and therefore, we 
consider that it is worth presenting them separately. 
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The paper is organized as follows. In Section [2] we prove the asymptotic estimates of W in L q 
norms for finite q and use them to study the behavior of the solution ul of the linear equation 
(jl.2p with a nonnegative, bounded initial datum satisfying (jl.5p with < a < N. 

In particular, we prove that the solution ui satisfies 
(1.9) t a / 2 \\uL{x, t) — ua{x, t)\\L°°(R N ) ~ * as t — )■ oo, 

where «a is the solution of the heat equation with diffusivity a and initial condition ul(x,0). 

Moreover, in this section we establish the precise behavior of ul for \x\ — > oo. 

In Section [3] we prove the global existence and uniqueness of the solution to (|1.1|) -with a more 
general absorption term — |u| p_1 ti- for our nonintegrable initial data with power like behavior at 
infinity. We also establish that the solution is bounded by B(l + |a;|) _a for a certain constant 
B. Then, when uq is nonnegative, we deduce that the solution is nonnegative. Finally, we prove 
that for every to > 

\x\ a u(x,to) — > A as \x\ — > oo. 

Finally, in Section U] we prove our main result. That is, we prove that (|1.6p holds when a < N 
and (|1.8p when a = N with u the solution to (jl.ip . 



2. Preliminary results and homogeneous equation 

In this section we study the asymptotic behavior of bounded solutions to the homogeneous 
equation 

u t = J(x-y)(u(y,t) -u(x,t))dy in R N x (0, oo), 

(2.1) J 

u(x,0)=u (x) mR N . 

Existence of bounded solutions for bounded initial data is a well know result that can be 
seen by different arguments. One possible construction is to convolve the initial datum with 
the fundamental solution (see (|2.2p below). Uniqueness of bounded solutions follows from the 
comparison principle for bounded solutions. (See [16] . Proposition 2.2 with 6 = for the case of 
smooth solutions. In order to get the result in our case we just have to approximate the initial 
datum by smooth functions that are uniformly bounded). 

In this section we first obtain estimates in L q (M. N ) for the "good" part W(x,t) of the funda- 
mental solution of the nonlocal diffusion equation. 

Then, we get the result on the asymptotic behavior as time goes to infinity of the solution of 
the equation for initial data that behave as a negative power at infinity. 

Finally, we prove that the solution behaves as the initial datum for every t > 0. 

Let us recall (see [S]) that the fundamental solution of the equation (|1.2p is given by 

(2.2) U(x,t) =e~ t 6 + W(x,t), 

where W is a smooth function. In fact, by taking Fourier transform in the space variables the 
authors show that W(-,t) belongs to the Schwartz class of C°° functions rapidly decaying at 
infinity together with their derivatives. 
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Moreover, in [13] by using the characterization of W in terms of its Fourier transform, the 
authors prove that 

(2.3) t N ' 2 \\W{x,t) - U a (x,t)\\ L o, iRN) < cr 1 ' 2 , 

where U a is the fundamental solution of the heat equation with diffusivity a. 

The Fourier transform method allows to get L q estimates only for q > 2. In order to get an 
L 1 estimate we observe that W is the solution to 

W t (x, t)= J J(x - y) (W(y, t) - W(x, t)) dy + e~ l J(x) 

W(x,0) = 0. 

In fact, let us call L the nonlocal operator. This is Lu = J * u — u. Then, 

Q 

W t (x, t) - LW(x, t) = U t (x, t) - LU(x, t) - T^e-*J + e _t L<5(x) 

= e- t 5 + e- t (J*5-5)(x) 
= e - *J(x). 

Moreover, since U(x,0) = 5 there holds that W(x,0) = 0. 

In particular, by the comparison principle, W > 0. On the other hand, since the unique 
bounded solution to the homogeneous nonlocal diffusion equation with initial datum uq = 1 is 
the function u(x,t) = 1, the representation of this solutions gives 

1 = u (x, t) = e"*n (x) + J W(x - y, t)u (y) dy = eT l + J W(y, t) dy. 

Thus, / W(y, t) dy = 1 - e -t < 1 for every t. 

Now, since j U a (y,t) dy = 1 for every t, interpolation with inequality (|2.3p gives the estimate 

N+l 1 1 

(2.5) \\W(x,t)-U a (x,t)\\ Lql(RN) <C q t ^ if + =i. 
In particular, for every 1 < q < oo, 

_N_ 

(2.6) \\W(x,t)\\ Lq , (RN) <C q t ^. 

With this estimate we can prove that for large t the solution u to (|2.1|) behaves as the solution 
of the heat equation with diffusivity a. In fact, 

Theorem 2.1. Let uq £ L°° , nonnegative such that there exist < a < N and A > with 
\x\ a uo(x) — >• A as \x\ — > oo. Let u be the solution of the equation ut = Lu with initial condition 
uq. Then, for every < [i < ^ there exists a constant such that 

(2.7) t a ' 2 \\u(x,t) - u A (x,t)\\ Lao{RN) < Cut-". 

Ln particular, if a < N, u(x, t) < Ct~ a / 2 and 
t a / 2 \\u(x,t)-U aA (x,t)\\ 
where U a) A is the solution to (|1.7|) . 
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Ifa = N, u{x,t) < Ct- N l 2 \ogt and 



t N/2 



u(x, t) 



logt 



C A ,NU a (x,t) 



L°°(\x\<K^i) 



as t — y oo \/K > 0, 



where U a is the fundamental solution of the heat operator with diffusivity a and Ca,n is a constant 
that depends only on A and N . 

Proof Since u 6 L*(R N ) if qa > N, by taking q = ^j^, e > 0, 

W(x - y, t)u (y) dy - J U a (x - y, t)u (y) dy 

<c q \\ Uo \\ q r^ = c e t-f c 1 -*)-^ 1 -*) 



r-i-t o yi-r N j 



C F t~2t~2N +S 2 



Cut-it-", 



with < fx < arbitrary. 
Thus, 

\\u(-,t) - u A (-,i)||oc < e~*||u |U + C M * _f * _M ' 

and (|277]) holds. 

The asymptotic behavior of u follows immediately from (|2,7p together with (|1.6p and (|1 .8j) 
respectively. □ 



Finally, we prove that u behaves at each positive time as its initial datum. 

Proposition 2.1. Let u be the solution of (|2.ip with an initial datum uq G L°° smc/i i/iai 
|x| a no(a;) — >■ ^4 > as |sc| — >■ oo /or some a > 0. Then, 

(1) For every T > i/iere exists a constant Ct > suc/i i/ioi < nrpraja */0 <t<T. 

(2) |x| Q u(x,i) — )• ^4 as | x| — >■ oo uniformly for t bounded. 

Proof. The proof follows from a fixed point argument. In fact, u is a fixed point of the operator 
Tv(x,t) = e^uoix) + y y e~^ s \j(x - y)v(y,s)dyds. 

Reciprocally, every bounded function that is a fixed point of the operator T in W N x (0,£o) 
is a bounded solution of problem (|2.ip . By uniqueness, it coincides with u in that time interval. 

Thus, in order to prove the proposition we will show that T has a fixed point in the set 

K := {v G L°°(R N x (0,t ))/(l + M) a |t;(x,i)| < 2B,\x\ a v{x,t) -»• A as |z| -> oo 
uniformly with respect to i G [0, to]}- 

The time step to will be independent of the constant B where B is a bound of (1 + |x|) a |ito(a;)|. 
Therefore, we can proceed inductively and find a constant for any time interval thus proving 
(1). Moreover, (2) follows immediately from this argument. 
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So, let v S K, and let w = Tv. Then, since the support of J is contained in a ball Br and 
\ x — y\ < R implies that \x\ < R + \y\, by assuming that R > 1 we get 

(l + |s|) a |iu(x,t)| < e- t (l + \x\) a \u (x)\ + £ J J(x - y)(l + \x\) a \v(y, s)\ dy ds 

< e- l B + J' J e~^J(x - y )(l + (R+ \y\) a )\v(y, s)\ dy ds 

< B + (R+ l) a Bt < 2B. 
If t < t = (R + l)~ a . 

Now, let us prove that \x\ a w(x,t) — > A as |x| — > oo uniformly for t G (0, to). In fact, if 
\x — y\ < R, there holds that \y\ > \x\ — R and 



1 



R_ R 

\y\ \x\ — R 



Thus, if \x\ is large enough, — l| < e. 
Therefore, 



< 



J(x - y)\x\ a v(y,s)dy - A 

< / J(x -y)\\y\ a v(y,s) - A 



dy + 



J(x - y)\\x\ a - \y\ a \ \v(y,s)\ dy 



< 



| W |>|x|-l 



J(x-y)edy + e / J(x - y)\y\ a \v(y, s)\ dy < (1 + 2B)e, 



if | a; | is large enough. 
Finally, 



\x\ a w{x,t) - A\ <e -t (|s| uo(a?) -A\ + f 

Jo 

<2(1 + B)e. 



-(*-«) 



J(x - y)\x\ a v(y,s)dy - A 



ds 



This ends the proof. 



□ 



3. Existence, uniqueness and first properties of the solution 

In this section we prove existence, uniqueness and time local properties of the solution to the 
nonlocal diffusion equation with absorption 

(3.1) ut = Lu — \u\ p ~ 1 u. 

In particular, we are interested in nonnegative bounded initial data that behave as a negative 
power at infinity. That is, which satisfy (jl.5p . 

Local existence and uniqueness of the solution to (|3.1|) with p > 1 and bounded initial datum 
follows by a fixed point argument. For instance, if we call S(t) the semigroup associated to 
the equation m = Lu in L°°(M> N ), we can find the solution as a fixed point in 1C := {v G 
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L°°(R N x (0,* ))/Woo <2||ti ||oo} of the operator 

Tv = S(t)uQ — / S(t — s)\v(x, s)\ p ~ l v(x, s) ds. 



t 

JO 

In fact, it is easy to see that for small to the operator T is a contraction from 1C to /C. 
Thus, there exists a bounded solution in the time interval (0, to). Since positive constants 
are supersolutions and negative constants are subsolutions to problem (|3.ip . and the function 
r |r| p_1 r is locally Lipschitz, the comparison principle (see, for instance, [16]) implies that 
the fixed point u is bounded by ||ito||tx>- Therefore, the solution can be extended for all times. 

Moreover, when uq > 0, the solution ul of the homogeneous equation (|1.2p with initial datum 
uq is nonnegative. Thus, ul is a supersolution and is a subsolution to (|3.1|) . By the comparison 
principle we deduce that 

(3.2) <u(x,t) <u L (x,t). 

For the type of initial data we are interested in, much more can be said. In fact, 



Theorem 3.1. Let uq € L°°(R ), uq > be such that \x\ a uo(x) — >■ A > as \x\ — >■ oo with 
a, A > 0. Let p > 1 and let u be the solution to (jl.ip . Then, 

(1) For every T < oo, t/iere exists a constant Ct such that 

u{x,t) <C(l + \x\)~ a fort<T. 

(2) If a < N , there exists a constant C such that u(x,t) < Ct~ a l 2 . If a = N , there exists a 
constant C such that u(x,t) < Ct~ N ^ 2 log t. 

(3) For every t > 0, |x| Q u(x,t) — >■ j4 as |x| — )• oo uniformly for t in bounded sets. 

Proof. (1) and (2) follow immediately from the estimate (|3.2p above and the results for ul 
(Proposition O (I))- 

In order to prove (3) we use the variations of constants formula 

u(x, t) = ul(x, t) — / S(t — s)u p (x, s) dx. 



We already know that the first term in the right hand side has the correct limit uniformly for 
t bounded (Proposition 12.1( 2)). Thus, we have to prove that the second term goes to faster 
than 

We think of u p (x,s) as a nonnegative, bounded, initial condition that satisfies that, for a 
constant B that is independent of s G [0,T], there holds that \x\ pa u p (x, s) < B. By Proposition 
12.11 we know that 

\x\ ap S{t - s)u p {x, s)<C T if < s < t < T, 
for a certain constant Ct- Thus, 

\x\ a S(t-s)u p {x,s) = Ix^-V^Sit- s)u p (x,s) < CtIxI"^" 1 ) < e, 

if \x\ is large. Therefore, 



x\ a / S(t- s)u p (x,s)dx 



< eT if < t < T and \x\ is large. 

□ 



NONLOCAL DIFFUSION WITH ABSORPTION 



!) 



4. Asymptotic behavior for the equation with absorption 



In this section we prove our main result. Namely, that in the supercritical case, the solution 
to with a bounded nonnegative initial datum uq satisfying (|1.5p has the same asymptotic 
behavior as the one of the homogeneous heat equation with diffusivity o. 

oo with 



Theorem 4.1. Let u > 0, u G L°°(R iy ) be such that \x\ a u (x) ->■ A > as \x 
< a < N. Let p > 1 + — . Let u be the solution to (jl.ip with u(x,0) = uq{x). 

Then, if a < N 

t a/2 \\u(x,t) - U a: A(x,t)\\ Loc ^ x ^ K ^ l - ) -> as t -)• oo VA > 0, 

where U a ^A is the solution to (|1.7p . 
7/a = JV, 



t AT/2 



logt 



CA,NU a (x,t) 



^0 as t ^ oo VA' > 

i°°(|a;|<ifV?) 

where U a is the fundamental solution of the heat operator with diffusivity a and Ca,n is a 
constant that depends only on A and N . 



Proof. By Theorem 13. II we know that (1 + \x\) a u(x,to) < C(to). Thus, by Theorem 12. II we have 
that the solution ul to 

v t = Lv in x (t ,oo) 
v(x,t ) = u(x,t ), 
satisfies 

t a/2 \\u L (x,t) - U a , A (x,t)\\ L ooQ x \< Kv r t) -> as t^oo VA>0 if a < N, 



(4.1) 



t N/2 



u L (x,t) 



logt 



C A ,NU(x,t) 



as t -)■ oo VA > if a = N. 



L™(\x\<Ky/i) 

Thus, the theorem will be proved if we show that for every e > there exists to such that 
t a / 2 \\u(x, t) — ul(x, t)\\ L oo/ K N) < e for t > 2t . 



In order to prove this result we need to estimate 



u(x, t) — ul(x, t) 



S(t — s)u p (x, s) ds. 



t 



We begin with the case a < N. 

Let us estimate the integrand S(t — s)u p (x, s). On one hand, since u(x, s) < Cs~ a ^ 2 , by the 
maximum principle for the solutions of Vt = Lv, there holds that 

(4.2) < S(t - s)u p (x, s) < Cs~ ap l 2 . 



On the other hand 
S(t- 



s)w [x,s 



e- {t - s) u p {x,s) + 



W{x - y,t - s)u p (y,s) dy 
< e -(t- S ) s -fp + ||^(.,t- s )|| 3 ,|K(-, S )|| g . 
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As p > 1 we can take q = Then, \\W{-,t- s)\\ q > < C(t - s) 2 « = C(t - s)~a. 
On the other hand, \\u p (x, s)\\ Lq( ^N^ = s)||^ P9mJVV Thus, since g = ^ > 1, 



< \\e- s u (-)\\ Pg + / W(x-y, S )u (y)dy 



<Ce- s + C[u ] q7OO \\W{;s)\\ r , 



with l + ^ = ^ + | and [no] gj00 = sup A>0 ^A 9 |{|uo| > A}|^ (see, for instance [TT], Theorem 
1.4.24). 

Since uq is in L q, °°, 



and hence, 



Therefore, 



|« p (-,s)||f < Ce" s + Cs-^ = Cs'V 1 , 



(4.3) 5(i - a) n p (x, s) < Ce"^^^ + C(t - s)~f s"f < C(t - a)" 



2 S 2 



Now we estimate the approximation error for t > 2tQ. We need to separate the integral and 
use estimate (|4.2|) for s G (t/2,t) and estimate (|4.3p for s G (to,t/2). We obtain, by using that 
§(p — 1) — 1 > (critical or supercritical cases), 

t /-t/2 ft 

S(t- s)u p (x,s)ds < C (t-s)"f S -f( p " 1 )ds + C / s-^ds 

to J t Jt/2 

*/2 „, /■* 



< (t/2)~f / s^Hds + c / s-fPds 

J tn Jt/2 



to Jt/2 

Thus, if to is large enough and t > 2to, we get, by using that we are in the supercritical case 
p > 1 + 2/a, 

if f S(t - s) u p (x, s) ds < Ct~^ {p ~ 1)+1 + CT § < e. 

Jt 



We remark that this is a sharp estimate that shows, in particular, that in the critical case the 
absorption and the diffusion terms are of the same order. 

Now we analyze the case a = N. 

We proceed as before. Estimate (|4.2p is changed to 

(4.4) < S(t - s)u p {x, s) < Cs- Np > 2 log p s. 



For the equivalent to estimate (|4.3p we have to proceed differently since the general form of 
Young's inequality that we have used is not valid when q = 1 as would be now the case. 
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We take instead q = with 5 > small to be chosen later (we could have proceeded in this 
way before but we wouldn't have gotten the sharp estimate), and use Young's inequality to get 

IKMIIi! = H;s)f LP < Ce~ ps + \\W(;s)\\ p Lr \\u Q \\ p Lq , 

with 1 + i = i + i. So that, £ = 1 - S - ± and g£ = f (p - 1 - 8p). Thus, 

\\u p (-,s)\\ L i < Ce~ ps + C^s-^p- 1 ^ < Cp8-% (P-V+», 
with ;U = y<5p as small as needed. 
Therefore, 

S(t - s)u p (x, s) < e- ( - t - s h-% p log p s+ \\W(-,t- s)|UIK0, 
(4.5) < C e - (i - s) s~f p logf s + C M (t - s)-T S -f (P-i)+*< 

< C/j( i_ s) -f s -f(p-D+M. 

As before, we use ()4.4p and (|4.5p to estimate the error. We have, by taking < \i < 
f(p-l)-l, 

||u(x,t) - UL(x,t)\\oo < || / 5(t - s)u p (x, s) ds\\oo 

ft/2 







<C M / {t-s)-^s-^ p - 1)+ ^ds + C f s-^ p \og p sds 

Jt Jt/2 

<c,{t/2y^t-^ p - 1)+1+ " + ct-TP+±i og p t. 

So that, if to is large enough and t > 2to we get 

t N / 2 \\u(x,t) - u L (x,t)\\oo < Cfo* +CfT (P-^logPt < e. 

So, the theorem is proved. □ 
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ASYMPTOTIC BEHAVIOR FOR A NONLOCAL DIFFUSION EQUATION 
WITH ABSORPTION AND NONINTEGRABLE INITIAL DATA. THE 

SUPERCRITICAL CASE 



JOANA TERRA AND NOEMI WOLANSKI 



Abstract. In this paper we study the asymptotic behavior as time goes to infinity of the 
solution to a nonlocal diffusion equation with absorption modeled by a powerlike reaction — u p , 
p > 1 and set in M. N . We consider a bounded, nonnegative initial datum uo that behaves like a 
negative power at infinity. That is, \x\ a uo(x) — > A > as \x\ — > oo with < a < N. We prove 
that, in the supercritical case p > 1 + 2/a, the solution behaves asymptotically as that of the 
heat equation -with diffusivity a related to the nonlocal operator- with the same initial datum. 



1. Introduction 

In this paper we study the asymptotic behavior as time goes to infinity of the solution to a 
nonlocal diffusion equation in R N with an absorption term. We are interested in the case in 
which the initial datum is nonintegrable. More precisely, we consider a bounded datum uq(x) 
such that |x| Q uo(3;) — > A > as |x| — > oo with < a < N. The problem under consideration is 
the following 



u t (x,t) = J J(x - y)(u(y,t) - u(x,t)) dy - u p (x,t) 
u(x, 0) = uq{x) 

where J € Cq°(^I n ), radially symmetric with J > and f J = 1. 

The term J J(x — y)(u(y,t) — u(x,t)) dy in (jl.ip represents diffusion. In fact, let u(x,t) be 
the density of a population at the point x at time t. The kernel J(x — y) may be seen as the 
probability distribution density of jumping from site y to site x. Thus, f J(x — y)u(y, t) dy 
represents the rate at which the population is arriving at the site x from all over space. By the 
symmetry of the kernel J, J J(x — y)u(x, t) dy represents then the rate at which it is leaving the 
point x. 

The absorption term —u p (x,t) in the equation represents a rate of consumption due to an 
internal reaction. 

This diffusion operator has been used to model several nonlocal diffusion processes in the last 
few years. See for instance [TJ [2J O HJ El US]- In particular, nonlocal diffusions are of interest 
in biological and biomedical problems. Recently, this kind of nonlocal operators have also been 
used for image enhancement |l(Jj . 



Key words and phrases. Nonlocal diffusion, boundary value problems. 
2000 Mathematics Subject Classification 35K57, 35B40. 
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It can be seen that -when properly rescaled- the nonlocal operator approximates aA, where 
A is the Laplace operator and a is a constant that depends only on J and the dimension N. 
(See, for instance, [8]). 

Another close relation to the heat operator with diffusivity o was found in [5J [13] where it 
was proven that, for bounded and integrable initial data, the asymptotic behavior as t tends to 
infinity of the solution ul to the equation without absorption 

(1.2) ut(x, t) = J J(x-y) (u(y, t) - u(x, t)) dy = Lu, 

is the same as that of the solution of the heat equation with diffusivity a and the same initial 
condition. Namely, 

(1.3) t N/2 \\u(x,t) ~ MU a (x,t)\\ LoaQxl < KVi) -»• ast^oo VK > 0, 

with U a the fundamental solution of the heat equation with diffusivity a and M = j uq the 
initial mass. The assumption on the initial condition is that uq G L°° n L . 

In a recent paper [T7j, the nonlocal equation with absorption (jl.ip was considered in the case 
u(x, 0) G L°° n L 1 and p > 1 + 2/N. The authors prove that n also satisfies (|1.3j) . Some results 
for p < 1 + 2/N can also be found in that paper. The critical case p = 1 + 2/N was left open. 

On the other hand, the asymptotic behavior of the solutions to the heat equation and of the 
heat equation with absorption are very well known for some nonintegrable initial data with a 
precise power like behavior at infinity. In fact, Kamin and Peletier in [T3] study the asymptotic 
behavior of the solution to 

u t = aAu-u p in R N x (0,oo), 

(1.4) 

u(x,0) =u (x) mR N , 
with no £ L°° and no > such that 

(1.5) |x| a no(x) — > A > as \x\ — > oo, 
with < a < N. In fact, the authors prove that 

(1.6) t a / 2 \\u(x,t)-U aA (x,t)\\ Lac(l , 

<KVt) ~ y as t y oo VK > 0, 

where U a ,A is the solution to 

U t = aAU - 5 P!a U p in R N x (0, oo), 

Here 5 PtOC = in the supercritical case p > 1 + 2/ a and 5 PtOC = 1 in the critical case p = 1 + 2/ a. 

Then, in [15] . Kamin and Ughi study the case where a = N in (|1.5|) and p > 1 + 2/N. They 
prove that 



-+0 as t -+ oo \/K > 0, 

L°°(\x\<K VI) 



(1.8) t^ 2 ^-C^^^t) 

logt 

with Ca,n a constant that depends only on A and N and J7 the fundamental solution of the 
heat operator with diffusivity a. 

Later on, in [12] . Herraiz fully analyzes problem (|1.4|) and presents a complete characterization 
including the asymptotic behavior outside the parabolas |x| < K\ft. 
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It is the purpose of this paper to prove that, for initial data satisfying (|1.5p with < a < N 
and p > 1 + 2/a, the asymptotic behavior of the solution to the nonlocal diffusion problem with 
absorption (jl.lj) is the same as the asymptotic behavior of solutions to the heat equation with 
absorption (jl.4p . That is, that (|1.6p holds in the case a < N and (jl.8p holds in the case a = N 
with u the solution of (jl.ip . 

The proofs of these results rely strongly on estimates for the fundamental solution of the 
equation without absorption (|1.2p . 

In [5] the authors observed that the fundamental solution of the operator in (jl.2p can be 
written as 

(1.9) U(x,t) =e- t S + W(x,t), 

where 5 is the Dirac measure and W(x,t) is a smooth function. 

Then, in [T3] the authors established L°° estimates for W and proved that 

t N/2 \\W(x,t) - U a (x,t)\\ La c(^ N) as t ->• oo, 
where U a is the fundamental solution of the heat equation with diffusivity a. 

This estimate allowed them to prove that the solution of the homogeneous nonlocal diffusion 
equation (|1.2p with bounded and integrable initial data behaves -for t tending to infinity- as the 
solution of the heat equation with diffusivity o, removing in this way the hypothesis that the 
Fourier transform of uq be in L 1 as was assumed in [5]. 

In the present paper, as we are considering nonintegrable initial data, these L°° estimates 
cannot be used. So, the first step is to get estimates in L 9 — norms for 1 < q < oo. 



As in [T7], the philosophy of the study of the asymptotic behavior for the problem with 
absorption in the supercritical case is that, as time tends to infinity, the absorption term goes 
to so fast that the behavior is the same as that of the equation without absorption. Thus, 
what we need to prove is essentially that, as time tends to infinity, the difference between the 
solutions to both problems goes to faster than each term separately. 

In order to prove this result it is necessary to compare the solution u of (jl.ip with the solution 
ul of the equation without absorption (jl.2p that coincides with u at a large time to- Once we 
prove that this difference is small the result follows if we know that the asymptotic behavior of 
ul is the same as that of f)1.4|) . To this end, we first analyze the asymptotic behavior for the 
nonlocal equation without absorption (jl.2p with an initial datum satisfying (|1.5p . In order to 
finish the argument, it is essential to know, and we prove it in this paper, that the behavior in 
space of the solution of our problem (jl.ip at any time to is the same as that of the initial datum 
Uq. That is, that u(-,to) satisfies (|1.5p at any positive time to. 

The techniques of this paper do not work in the critical case p = 1 + 2/a. In fact, for such p, 
the asymptotic behavior is that of the heat equation with absorption since the diffusion and the 
absorption terms are of the same order. We have to proceed in a completely different way by 
going back to the techniques of the original paper for the heat equation p~3] , and we do so in [18] . 
The method in [18] also works in the supercritical case considered in this paper. Nevertheless, 
the ideas in the present paper are simpler and rather easy to implement and therefore, we 
consider that it is worth presenting them separately. 
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The paper is organized as follows. In Section [2] we prove the asymptotic estimates of W in L q 
norms for finite q and use them to study the behavior of the solution ul of the linear equation 
(jl.2p with a nonnegative, bounded initial datum satisfying (jl.5p with < a < N. 

In particular, we prove that the solution ui satisfies 

(1.10) t a ' 2 \\uL(x, t) — UA(x,t)\\ Lao mN) — > as t — y oo, 

where «a is the solution of the heat equation with diffusivity a and initial condition ul(x,0). 

Moreover, in this section we establish the precise behavior of ul for \x\ — > oo. 

In Section [3] we prove the global existence and uniqueness of the solution to (jl.ip -with 
a more general absorption term —\u\ p ~ 1 u- for our nonintegrable initial data with power like 
behavior at infinity. We also establish that the solution is bounded, in any time interval [0,T], 
by Ct(1 + |^|) _a for a certain constant Ct- Then, when uq is nonnegative, we deduce that the 
solution is nonnegative. Finally, we prove that for every to > 

|x| a u(x,to) A as \x\ — > oo. 

Finally, in Section U] we prove our main result. That is, we prove that (|1.6|) holds when a < N 
and (|1.8|) when a = N with u the solution to (|1.1|) . 



2. Preliminary results and homogeneous equation 

In this section we study the asymptotic behavior of bounded solutions to the homogeneous 
equation 

u t = J(x -y)(u(y,t) -u(x,t)) dy in R N x (0, oo), 

(2.1) J 

u(x,0)=u (x) mR N . 

Existence of bounded solutions for bounded initial data is a well know result that can be 
seen by different arguments. One possible construction is to convolve the initial datum with 
the fundamental solution (\1.9\i . Uniqueness of bounded solutions follows from the comparison 
principle for bounded solutions. (See [16] . Proposition 2.2 with 6 = for the case of smooth 
solutions. In order to get the result in our case we just have to approximate the initial datum 
by smooth functions that are uniformly bounded). 

In this section we first obtain estimates in L q (W N ) for the "good" part W(x,t) of the funda- 
mental solution of the nonlocal diffusion equation. 

Then, we get the result on the asymptotic behavior as time goes to infinity of the solution of 
the equation for initial data that behave as a negative power at infinity. 

Finally, we prove that the solution behaves as the initial datum for every t > 0. 

Let us recall (see [5]) that the fundamental solution of the equation (jl.2p is given by 

(2.2) U(x,t) =e~ t 8 + W(x,t), 

where W is a smooth function. In fact, by taking Fourier transform in the space variables the 
authors show that W(-,t) belongs to the Schwartz class of C°° functions rapidly decaying at 
infinity together with their derivatives. 
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Moreover, in [13] by using the characterization of W in terms of its Fourier transform, the 
authors prove that 

(2.3) t N l 2 \\W{x,t) - U a (x,t)\\v*>&») < Ct~ 1/2 , 

where U a is the fundamental solution of the heat equation with diffusivity a. 

The Fourier transform method allows to get L q estimates only for q > 2. In order to get an 
L 1 estimate we observe that W is the solution to 

W t (x, t) = J J(x - y) (W(y, t) - W(x, t)) dy + e~* J{x) 

W(x,0) = 0. 

In fact, let us call L the nonlocal operator. This is Lu = J * u — u. Then, 

W t (x, t) - LW(x, t) = UAx, t) - LU(x, t) - 4-e~*<5 + e^L5(x) 

at 

= e- t S + e- t (J*6-S)(x) 
= e-*J(x). 

Moreover, since U(x, 0) = 8 there holds that W(x,0) = 0. 

In particular, by the comparison principle, W > 0. On the other hand, since the unique 
bounded solution to the homogeneous nonlocal diffusion equation with initial datum uq = 1 is 
the function u(x,t) = 1, the representation of this solutions gives 

l = u(x,t) = e~ t u (x)+ / W(x-y,t)u (y)dy = e~ t + / W(y,t)dy. 



Thus, / W(y, t) dy = 1 - e~* < 1 for every t. 

Now, since j U a (y,t) dy = 1 for every t, interpolation with inequality (|2.3|) gives the estimate 

N+l 1 1 

(2.5) \\W(x,t)-U a (x,t)\\ Lq/iViN) <C q t ^ if - + - = 1. 

In particular, for every 1 < q < oo, 

(2.6) \\W(x,t)\\ Lq , (RN) <c q r%. 



With this estimate we can prove that for large t the solution u to (|2.ip behaves as the solution 
of the heat equation with diffusivity a. In fact, 

Theorem 2.1. Let uq G L°°, nonnegative such that there exist < a < N and A > with 
\x\ a uo(x) — >• A as \x\ — > oo. Let u be the solution of the equation ut = Lu with initial condition 
uq. Then, for every < fi < ^ there exists a constant such that 

(2.7) t a ' 2 \\u(x,t) - u A (x,t)\\ Laa{RN) < Cpt-*. 

Ln particular, if a < N, due to the behavior of ua (see, for instance [H]J there holds that 
u(x,t) < Cr a / 2 and 

t^ 2 \\ U (x,t)-U aA (x,t)\\ Lao(M 
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where U aj A is the solution to (|1.7p . (Note that the results of [141 115] apply also to the solution 
ua of the homogeneous heat equation. In this case, the constant 5 ay A in (jl-7j) is zero so that, 
U a ,A is a solution to the homogeneous heat equation). 

If a = N, due to the results in [15] for the solution of the homogeneous heat equation, there 
holds that u(x,t) < Ct~ N ^ 2 logt and 

t N/2 ™^)__ c AN UJx,t) ^0 ast^oo VK>0, 

where U a is the fundamental solution of the heat operator with diffusivity a and Ca,n is a constant 
that depends only on A and N . 

Proof. Since uq € L q (M. N ) if qa > N, by taking q = ^jr^, e > and applying (|2.5|) we get, 

W(x - y,t)u Q (y)dy - / U a (x - y,t)u (y) dy 



< c q \\u \\ q r^ = c e rf d-^-wd-) 
= c B t-St-^ +e 2( 1+ *) =c^r\ 

with < fx < ^ arbitrary. 
Thus, 

IK-,*) -u A (-,t)\\L°°(BL N ) < e_ *ll u o||L-(RJV) + C l jr^t~ lx , 

and (O) holds. 



The asymptotic behavior of u follows immediately from (|2.7|) together with (|1.6|) and (|1.8|) 
respectively applied to the solution ua of the homogeneous heat equation (see [14] [15] ) . □ 

Finally, we prove that u behaves at each positive time as its initial datum. 

Proposition 2.1. Let u be the solution of ()2.1|) with an initial datum uq € L°° such that 
\x\ a uo(x) — > A > as | x | — >■ oo for some a > 0. Then, 

(1) For every T > there exists a constant Ct > such that \u(x,t)\ < mq^na */0 — * — 

(2) |x| Q w(x,i) —)• j4 as | x| oo uniformly for t bounded. 

Proof. The proof follows from a fixed point argument. In fact, u is a fixed point of the operator 
Tv(x,t) = e _t uo(x) + / / e~^~ s \j(x-y)v(y,s)dyds. 



Jo J 

Reciprocally, every bounded function that is a fixed point of the operator T in W N x (0,to) 
is a bounded solution of problem f)2. 1 [) . By uniqueness, it coincides with u in that time interval. 

Thus, in order to prove the proposition we will show that T has a fixed point in the set 

K := {v E L°°(R N x (0,t ))/(l + \x\) a \v{x,t)\ < 2B,\x\ a v{x,t) -> ,4 as |z| -> oo 

uniformly with respect to i S [0, to]}- 

The time step to will be independent of the constant B where B is a bound of (1 + |x|) a |ito(a;)|. 
Therefore, we can proceed inductively and find a constant for any time interval thus proving 
(1). Moreover, (2) follows immediately from this argument. 
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So, let v S K, and let w = Tv. Then, since the support of J is contained in a ball Br and 
\ x — y\ < R implies that \x\ < R + \y\, by assuming that R > 1 we get 

(1 + \x\) a \w(x,t)\ < e"*(l + |x|) a |uo(a:)| + J J e~^J(x - y)(l + \x\) a \v(y, s)\ dy ds 

< e^B + J* J e-^J(x -y)(l + (R+ \y\) a )Hy, s)\ dy ds 

< B + (R + l) a Bt < 2B. 
If t < t = (R + l)~ a . 

Now, let us prove that \x\ a w(x,t) — > A as |x| — > oo uniformly for t £ (0, to). In fact, if 
\x — y\ < R, there holds that \y\ > \x\ — R and 



fl 1 

\y\ 



R R 

< M < 1 i 5" 

\y\ \x\ — R 



Thus, if \x\ is large enough, — l| < e. 
Therefore, 

J(x - y)\x\ a v(y,s)dy - A 
< [ J(x-y)\\y\ a v 
< 



< 



dy + 



J(x - y)\\x\ a - \y\ a \ \v(y,s)\ dy 



[ J(x-y)edy + e [ J(x - y)\y\ a \v(y, s)\ dy < (1 + 2B)e, 

J\y\>\x\-R J 



if | a; | is large enough. 
Finally, 



\x\ a w{x,t) - A\ <e^(\x\ a u {x) - A ) + / 

Jo 

<2{l + B)e. 



J(x - y)\x\ a v(y,s)dy - A 



ds 



This ends the proof. 



□ 



3. Existence, uniqueness and first properties of the solution 

In this section we prove existence, uniqueness and time local properties of the solution to the 
nonlocal diffusion equation with absorption 

lit = Lu — \u\ p ~ 1 u 
(3.1) 1 ' 

u(x, 0) = uq(x) 

In particular, we are interested in nonnegative bounded initial data that behave as a negative 
power at infinity. That is, which satisfy (|1.5|h 

Local existence and uniqueness of the solution to (|3.ip with p > 1 and bounded initial datum 
follows by a fixed point argument. For instance, if we call S(t) the semigroup associated to 
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the equation ut = Lu in L°°(fil N ), we can find the solution as a fixed point in K, := {v £ 
L°°(R N x (0, t )) / Halloo < 2||ii ||oo} of the operator 

ft 



Tv = S(t)uQ— / S(t — s)\v(x, s)\ p 1 v(x,s)ds. 
Jo 



In fact, it is easy to see that for small to the operator T is a contraction from KL to JC. 
Thus, there exists a bounded solution in the time interval (0, to). Since positive constants 
are supersolutions and negative constants are subsolutions to problem (|3.ip . and the function 
r — > \t\ p ~ 1 t is locally Lipschitz, the comparison principle (see, for instance, [IB]) implies that 
the fixed point u is bounded by ||«o||oo- Therefore, the solution can be extended for all times. 

Moreover, when uq > 0, the solution ul of the homogeneous equation (|1.2p with initial datum 
ur, is nonnegative. Thus, ul is a supersolution and is a subsolution to (|3.ip . By the comparison 
principle we deduce that 

(3.2) < u(x,t) < u L (x,t). 

For the type of initial data we are interested in, much more can be said. In fact, 

Theorem 3.1. Let uq G L od (M. n ), no > be such that \x\ a UQ(x) — >■ A > as \x\ — )• oo with 
a, A > 0. Let p > 1 and /et tt be the solution to (jl.ip . Then, 

(1) For every T < oo, t/iere exists a constant Ct such that 

u{x,t) <C T (l + \x\)- a fort<T. 

(2) If a < N , there exists a constant C such that u(x,t) < Ct~ a l 2 . If a = N , there exists a 
constant C such that u(x,t) < Ct~ N l 2 log t. 

(3) For every t > 0, |x| Q n(x,t) — >• yl as |x| — ?■ oo uniformly for t in bounded sets. 

Proof. (1) and (2) follow immediately from the estimate (|3.2p above and the results for ul 
(Proposition O(l)). 

In order to prove (3) we use the variations of constants formula 

u(x,t) = UL(x,t) — / S(t — s)u p (x, s) ds. 



We already know that the first term in the right hand side has the correct limit uniformly for 
t bounded (Proposition 12.1( 2)). Thus, we have to prove that the second term goes to faster 
than 

We think of u p (x,s) as a nonnegative, bounded, initial condition that satisfies that, for a 
constant B that is independent of s £ [0,T], there holds that \x\ pa u p (x, s) < B. By Proposition 
12.11 and Theorem 12.11 we know that 



\x\ ap S(t - s)u p (x, s)<C T if < s < t < T, 
for a certain constant Ct- Thus, 

\x\ a S(t- s)u p {x,s) = \x\- a( P-V\x\ ap S(t- s) U P(x,s) < C T \x\- a( > p -V <£, 
if |x| is large. Therefore, 



S(t — s)u p (x, s) dx 



o 



< eT if < s < t < T and Ixl is large. 
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4. Asymptotic behavior for the equation with absorption 



In this section we prove our main result. Namely, that in the supercritical case, the solution 
to with a bounded nonnegative initial datum uq satisfying (|1.5|) has the same asymptotic 
behavior as the one of the homogeneous heat equation with diffusivity a. 

Theorem 4.1. Let uq > 0, uq 6 L°°(R ) be such that \x\ a uo{x) — > A > as \x\ — > oo with 
< a < N. Let p > 1 + ^. Lei u 6e £/ie solution to (jl.ip to£/j u(a;, 0) = uo(x). 

Then, if a < N 

t a l 2 \\u{x,t) - U a , A (x,t)\\ 

where U a ^ A is the solution to (jl.7p , 
//a = JV, 

-j 7 C A ,NUa(X,t) 

log£ ' L°°(|z|<i<Vt) 

where U a is the fundamental solution of the heat operator with diffusivity a and Ca,n is a 
constant that depends only on A and N . 



t N/2 



as t ->■ oo MK > 



Proof. By Theorem 13. II we know that (1 + \x\) a u{x,to) < C(£o). Thus, by Theorem 12. II we have 
that the solution ul to 

v t = Lv in R N x (£ , oo) 
v(x,t ) = u(x,t ), 
satisfies 



(4.1) 



t a l 2 \\u L {x,t) - U a>A {x,t) 



\L°°{\x\<Ky/i) 



t N/2 



u L (x,t) 



log* 



-)■ as £ -> oo VET > if a < iV, 



^ as £ ^ oo Vif > if a = N. 



L°o(\x\<KVi) 

Thus, the theorem will be proved if we show that for every e > there exists £o such that 
t a l 2 \\u{x, t) — Ul(x, £)||Loo( ffi jv) < e for £ > 2£o- 
Here ul is the solution of (|4.ip . 

In order to prove this result we need to estimate 



u(x, t) — ul(x, t) 



S(t — s)u p (x, s) ds. 



t 



We begin with the case a < N. 

Let us estimate the integrand S(t — s)u p (x, s). On one hand, since by (|3.2p and Theorem 12.1 
there holds that u(x, s) < Cs~ a / 2 , the maximum principle applied to the solutions of vt = Lv, 
renders the estimate 

(4.2) < S{t - s)u p {x, s) < Cs~ ap/2 . 



10 



J. TERRA AND N. WOLANSKI 



On the other hand, 

S (t - .) .) = ,-<->„<*, •) + / W (, - y, t - .W, .) * 
< e -»-«)»-f +||H^(-,(-»)||^||«f(-,,)||,. 

_N_ a 

As j> > 1 we can take g = f . Then, ||W(-,t - s)|| g / < C(i - s) 2 « = C(i - s)" . 
On the other hand, s)||l<j(irjv) = \\u(x, s )\\ P LP q(^Ny Thus, since q = ^ > 1, 

\\u p (;s)\\l < ||e- s u (-)IU+ / W(x-y,s)u (y)dy 

J pq 

<Ce- s + C[u ] q}OO \\W{;s)\\ r , 

with l + ^ = ^ + | an d [wo]g,oo = sup A>0 ^A 9 |{|no| > A}|^ (see, for instance [TT], Theorem 
1.4.24). 

Since uq is bounded and |x| a tio(x) — > A as \x\ — > oo there holds that uq G L q '°° and, 



and hence, 



\u p (-,s)\\l < Ce~ s + Cs~^ = Cs~^ }1 , 



Therefore, 

(4.3) S(t-s)u p (x,s) < Ce-^s-^P + C(t - s)-%s-^ p -V < C(t - a)~f a~f 

Now we estimate the approximation error for i > 2io- We need to separate the integral and 
use estimate (|4.2|) for s G (t/2,t) and estimate (|4.3p for s G (to,t/2). We obtain, by using that 
§(p — 1) — 1 > (critical or supercritical cases), 

/ S(i-s)u p (a;,s)ds < C/ {t - s)"^ s~^ p - l) ds + C s'^ds 

J t Jt Jt/2 

ft/2 ft 

< (t/2)"t / s-fCP-^ds + C / s~f p ds 

A) Jt/2 

< ct-ft - f(p - 1)+1 + c^fp +1 . 

Thus, if to is large enough and t > 2to, we get, by using that we are in the supercritical case 
p > 1 + 2/a, 

tf /* s(* - S ) ^(x, s ) d s < cto f (p " 1)+1 + cr f < e. 

We remark that this is a sharp estimate that shows, in particular, that in the critical case the 
absorption and the diffusion terms are of the same order. 



Now we analyze the case a = N. 
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We proceed as before. Estimate (|4.2p is changed to 
(4.4) < S(t - s)u p (x, s) < Cs~ Np/2 log p s. 



For the equivalent to estimate (|4.3p we have to proceed differently since the general form of 
Young's inequality that we have used is not valid when q = 1 as would be now the case. 

We take instead q = with 5 > small to be chosen later (we could have proceeded in this 
way before but we wouldn't have gotten the sharp estimate), and use Young's inequality to get 

KMIUi = IK,* < Ce~ ps + \\W(;s)\\l r \\u Q \\l q , 
with 1 + i = i + i. So that, £ = 1 - 6 - \ and = f (p - 1 - 8p). Thus, 
K(-,s)|| L i < Ce~ ps + C^s-tCp- 1 ^ < C^a-T (p-i)+f, 

with fi = y$P as small as needed. 
Therefore, 

S(i - S)U P (X,S) < e- ( - t - S h-fn g P S+ \\W(;t- S)\U\U P (; S)^ 

(4.5) < Ce- (i - s) s-f p bgP s + C M (t - s )~f <TT (p-^+z* 

AT iV 



w w -(p-i)+/i 



As before, we use (|4.4p and (|4.5|) to estimate the error. We have, by taking < fj, < 
f(P"l)-l, 



\u(x,t) - ^(x.^Hioo^JV) < || / S 1 ^ - s)'U P (x,s)(is|| i cx ) ( M iV) 

ft/2 



<CW (t-s)-% S -% ip - 1)+ »ds + C s-% p log p sds 

J t Jt/2 



'to Jt/2 

N _JV (p _ 1) + 1+M 



< C^(t/2)~ T *o T +C*" fp+1 log p *. 
So that, since p > 1 + 2/N, if to is large enough and t > 2io we get 

So, the theorem is proved. □ 
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